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ABSTRACT 

In this manuscript, the structure of broad emission hne regions (BLRs) of weU- mapping 
double-peaked emitter (AGN with broad double-peaked low-ionization emission lines) 
3C390.3 is studied. Besides the best fitted results for double-peaked broad optical 
balmer lines of 3C390.3 by theoretical disk model, we try to find another way to 
further confirm the origination of double-peaked line from accretion disk. Based on 
the long-period observed spectra in optical band around 1995 collected from AGN 
WATCH project, the theoretical disk parameters of disk-like BLRs supposed by ellip- 
tical accretion disk model (Eracleous et al. 1995) have been well determined. Through 
the theoretical disk-like BLRs, characters of observed light-curves of broad double- 
peaked Ha of 3C390.3 can be well reproduced based on the reverberation mapping 
technique. Thus the accretion disk model is preferred as one better model for BLRs of 
3C390.3. Furthermore, we can find that different disk parameters should lead to some 
different results about size of BLRs of 3C390.3 from the one measured through ob- 
servational data, which indicates the measured disk parameters are significantly valid 
for 3C390.3. After that, the precession of theoretical elliptical disk-like BLRs being 
considered, we can find that the expected line profile in 2000 by theoretical model is 
consistent with the observed line profile by HST around 2000. Based on the results, 
we can further believe that the origination of broad double-peaked balmer emission 
lines of 3C390.3 are from accretion disk around central black hole. 
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1 INTRODUCTION 

The properties of objects with broad double-peaked low- 
ionization emission lines (hereafter, double-peaked emitters) 
have been studied for more than two decades, since the 
double-peaked emitters were found in nearby radio galaxies 
(Stauffer et al. 1983, Oke 1987, Perez et al. 1988, Halpern 
1990, Chen et al. 1989, Chen & Halpern 1989). Some the- 
oretical models have been proposed to explain the prop- 
erties of double-peaked broad emission lines which can be 
used as one probe of the broad-line regions of active galac- 
tic nuclei (Eracleous et al. 2009), such as binary black hole 
model (Begelman et al. 1980, Gaskell 1983, 1996, Boroson 
& Lauer 2009, Lauer & Boroson 2009, Zhang et al. 2007), 
double stream model (Zheng et al. 1990, 1991, Vellieux & 
Zheng 1991), accretion disk model (Chen et al. 1989, Chen 
& Halpern 1989, Eracleous et al. 1995, 1997, Bachev 1999, 
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HartnoU & Blackman 2000, 2002, Karas et al. 2001, Gezari 
et al. 2007, Flohic & Eracleous 2008, Lewis et al. 2010, Tran 
2010, Chornock et al. 2010, GaskeU 2010) etc.. Ahhough 
which theoretical model is more preferred for double-peaked 
emitters is still an open question, the accretion disk model 
which indicates double-peaked broad emission lines are from 
central accretion disk is so far the more widely accepted 
model applied to explain the properties of double-peaked 
broad emission lines. 

Due to the un-obscured emitting region on the receding 
jet, the double-stream model can be ruled out for double- 
peaked emitter 3C390.3 (Livio & Xu 1997). Through the 
long-period observational results, binary black hole model 
has also been ruled out due to unreasonable large cen- 
tral BH masses for some double-peaked emitters, such as 
Arpl02B, 3C390.3 (Eracleous et al. 1997). Furthermore, 
based on the long-period variations of double-peaked broad 
emission lines, accretion disk model is successfully applied 
to reproduce the characters of varying observed line prohles 
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of double-peaked emitters, such as the individual double- 
peaked emitter NGC1097 (Storchi-Bergmann et al. 1995, 

1997, 2003), and one sample of double-peaked emitters 
(Lewis et al. 2010, Gezari et al. 2007). However, we should 
note that although accretion disk model is so far the widely 
accepted theoretical model for double-peaked emitters, the 
other theoretical models can be expected to be preferred for 
some special double-peaked emitters, such as the x-shaped 
radio objects which could be candidates for objects with bi- 
nary black holes in central regions (Merritt & Ekers 2002, 
Zhang et al. 2007). Thus, only based on the double-peaked 
appearance of broad emission lines, there is no confirmed 
way to affirm which theoretical model is more preferred for 
double-peaked emitter. 

Among the sample of double-peaked emitters (Era- 
cleous & Halpern 1994, 2003, Eracleous et al. 1995, Strat- 
eva et al. 2003), 3C390.3 is one well studied double-peaked 
emitter, also one well-studied mapping AGN (Dietrich et al. 

1998, O'Brien et al. 1998, Leighly et al. 1997, Shapovalova 
et al. 2001, 2010, Sergeev et al. 2002, Peterson et al. 2004, 
Pronik & Sergeev 2007, Sambruna et al. 2009, Jovanovic 
et al. 2010). Based on the reverberation mapping technique 
(Blandford & Mckee 1982) and virialization method (Peter- 
son et al. 2004, Collin et al. 2006, Peterson & Bentz 2006, 
Peterson 2010), size of BLRs (distance between central black 
hole and broad line emission gas clouds) and central virial 
BH masses of 3C390.3 have been well determined (Peterson 
et al. 2004, Onken et al. 2004, Kaspi et al. 2005, Bentz et al. 
2006, 2009, Brandon & Bechtold 2007). Besides the common 
BH masses and size of BLRs, through the theoretical accre- 
tion disk model applied for observed double-peaked broad 
emission lines, some basic disk parameters of supposed the- 
oretical disk-like BLRs of 3C390.3 have also been measured 
(Eracleous & Halpern 1994), such as the inclination angle of 
accretion disk, the inner and outer radius of BLRs in accre- 
tion disk etc.. However, only due to the best fitted results for 
double-peaked balmer lines by accretion disk model, we can 
not firmly confirm that the double-peaked broad emission 
lines are exactly coming from disk-like BLRs in accretion 
disk around central black hole for 3C390.3. 

Besides the best fitted results for double-peaked broad 
emission lines by theoretical models, based on the pioneer 
work by Blandford & Mckee (1982), geometrical structures 
of BLRs can be mathematically structured by so-called 
transfer function in the reverberation mapping technique 
through some special mathematical methods (such as the 
Maximum Entropy Method, Narayan & Nityananda, 1986) 
applied to the properties of long-period observational vary- 
ing both continuum emission and broad line emission (Pe- 
terson et al. 1993, 1994, Home et al. 1991, Goad et al. 1993, 
Wanders & Home 1994, Pijpers & Wanders 1994, Krolik 
1994, Winge et al. 1995, Bentz et al. 2010). However, it is 
unfortunate that there should be NOT unique solution to 
the so-called transfer function used to determine the struc- 
tures of BLRs of AGN, due to the less complete informa- 
tion about variations of continuum and broad emission lines 
(Maoz 1996). Thus, even based on the solution of transfer 
function, the supposed theoretical and mathematical geo- 
metrical structures of BLRs of AGN can NOT still be firmly 
confirmed. 

In the previous remarkable work about structures of 
BLRs of AGN, the main objective is to determine the struc- 



tures of BLRs through observational properties/characters. 
However, we think it will be also very interesting to check 
whether the observational properties of long-period vari- 
ations of broad line emission and/or continuum emission 
could be better reproduced through the reverberation map- 
ping technique, based on the expected structures of BLRs 
supposed by theoretical model, which is the main objective 
of our paper. 

The structures of our paper are as follows. Section 2 
gives some information about observed spectra of 3C390.3 
around 1995 and how to find theoretical disk parameters for 
the supposed disk-like BLRs through the best fitted results 
for broad double-peaked balmer lines by theoretical ellipti- 
cal accretion disk model (Eracleous et al. 1995). Section 3 
shows the results through the reverberation mapping tech- 
nique under the disk-like BLRs for 3C390.3. Finally section 
4 gives some detailed discussions and conclusion. The cos- 
mological parameters Hq = 70km ■ s~^Mpc~^, f^A = 0.7 and 
f2m = 0.3 have been adopted here. 



2 FITTED RESULTS FOR SPECTRA OF 
3C390.3 AROUND 1995 

As one well mapping double-peaked broad 
line AGN, 3C390.3 (z=0.056) is one tar- 
g et included in the project of AGN W ATCH 
(httpi/Zwww.astronomy.ohio-state.edu/'agnwatch/l, which 
is a consortium of astronomers who have studied the inner 
structure of AGN through continuum and emission-line 
variability. From the AGN WATCH project, we can collect 
133 spectra of 3C390.3 in optical band observed around 
1995 by different instruments in different observatories. 
Meanwhile the light curves of continuum and broad emission 
lines after corrections of some necessary contamination are 
also collected from the website. The detailed descriptions 
about the instruments and techniques for the observed 
spectra can be found in Dietrich et al. (1998). Here, we do 
not describe the information any more. We mainly focus 
on how to determine the disk parameters by theoretical 
accretion disk model for 3C390.3. 

In order to obtain reliable disk parameters for disk-like 
BLRs of 3C390.3 supposed by the theoretical accretion disk 
model applied to fit double-peaked broad emission lines, we 
prefer to use double-peaked broad Ha rather than double- 
peaked broad H/?, because of the more apparent double 
peaks of broad Ha (one peak of double-peaked broad H/3 
is mixed by [OIII]A4959, 5007A doublet) and more stronger 
broad Ha than broad H/3. Thus we mainly focus on the 
66 out of the 133 observed spectra with reliable broad Ha 
within rest wavelength from 6200A to 7300A. 

There are so far several kinds of accretion disk mod- 
els which can be applied for double-peaked emitters, cir- 
cular with/without spiral arms accretion disk model(Chen 
et al. 1989, Chen & Halpern 1989, HartnoU & Blackman 
2002), elliptical accretion disk model (Eracleous et al. 1995), 
warped accretion disk model (HartnoU & Blackman 2000), 
stochastically perturbed accretion disk model (Flohic & Er- 
acleous 2008) etc.. In this paper, the elliptical accretion 
disk model (Eracleous et al. 1995) is preferred, because the 
model can explain most of the observational spectral fea- 
tures of broad double-peaked emission lines (especially fea- 
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tures for extended asymmetric line wings) with less num- 
ber of necessary model parameters. Furthermore, the most 
part of flux density of broad double-peaked line emission 
is from disk-like BLRs into accretion disk, the existence of 
arms (HartnoU & Blackman 2002) and/or warped structures 
(Hartnoll & Blackman 2000) and/or bright spots (Flohic & 
Eracleous 2008) are mainly applied for subtle structures of 
double-peaked line profiles (such as some cusps around the 
peaks etc.), which have few effects on the results based on 
Cross Correlation Function in reverberation mapping tech- 
nique. The detailed description of the elliptical accretion 
disk model with seven model parameters can be found in 
Eracleous et al. (1995). The seven necessary parameters are 
inner radius ro, out radius r\, eccentricity of elliptical rings 
e, local broadening velocity a, inclination angle of disk- 
like BLRs i, line emissivity slope (/r oc r"') and orienta- 
tion angle of elliptical rings 00 ■ Through the Levenberg- 
Marquardt least-squares minimization method applied to 
fit the double-peaked broad Ha by elliptical accretion disk 
model as what we have done for the x-shaped radio source 
SDSS J1130-I-0058 (Zhang et al. 2007), the seven disk pa- 
rameters in elliptical accretion disk model can be well de- 
termined. The procedures to determine the disk parameters 
for BLRs of 3C390.3 are as follows. 

First and foremost, in order to check the results from 
reverberation mapping technique under the structures of 
BLRs supposed by the accretion disk model, we should 
measure reliable theoretical disk parameters for disk-like 
BLRs of 3C390.3. Before starting to fit the observed broad 
double-peaked Ha, we first check the spectra of 3C390.3 
by eye, and find that the spectra marked with 'ce' (cor- 
responding observatory can be found in Table 1 in Di- 
etrich et al. 1998) have some unexpected absorption fea- 
tures around 6500A in rest wavelength. The weird absorp- 
tion features can not be found in the other spectra ob- 
served by instruments in other observatories, which should 
bring some larger uncertainties in measured disk param- 
eters when broad Ha is fitted by accretion disk model. 
Thus the absorption features (rest wavelength from 6507A 
to 6525A) in the 24 spectra marked with 'ce' should be re- 
jected, when the broad Ha are fitted by the elliptical accre- 
tion disk model for 3C390.3. Besides the weird absorption 
features in the 24 spectra marked with 'ce', in order to re- 
duce the effects of narrow emission lines ([OI]A6300, 6363A, 
[Nn]A6548,6583A, [Sn]A6716, 673lA and narrow Ha) in aU 
the spectra, the parts of narrow emission lines are also 
masked when broad Ha are fitted, rest wavelength from 
62681 to 6326 A for [OI]A6300A, rest wavelength from 6350A 
to 6388A for [OI]A6363A, rest wavelength from 6533A to 
6619A for [Nn]A6548, 6583A and narrow Ha, rest wave- 
length from 6703 A to 6743 A for [Sn]A6716, 6731 A (the grey 
shadows shown in Figure [Jl Figure [2] and Figure (5]). In other 
words, only the double-peaked broad component of Ha with- 
out weird absorption features and without narrow emission 
lines are fitted. After the prepared work, the double-peaked 
broad Ha are fitted two times as follows. 

Firstly, the 39 spectra with high quality (more than 
280 reliable data points within rest wavelength range from 
6I6OA to 6920A) and without unexpected absorption fea- 
tures around 6500A are fitted by the elliptical accretion disk 
model. Then through the measured disk parameters, the two 
parameters of eccentricity and inclination angle are firstly 



determined for the 39 spectra, because eccentricity and in- 
clination angle should not be changed with passage of time. 
The mean values of the two parameters for the 39 spectra 
are accepted as the reliable values for the two parameter. 
Then all the 66 spectra are fitted again by the accretion 
disk model with the accepted constant inclination angle and 
eccentricity. Figure [T] shows the best fitted results for some 
selected broad Ha by elliptical accretion disk model with 
MJD from 49770 to 50051. In the figure, two spectra with 
unexpected absorption features (marked with '49860ce' and 
'49984ce') are also shown, the shadow areas represent the 
masked ranges for narrow emission lines. Figure [2] shows the 
corresponding residuals {yobs — Vmodei , observed data minus 
expected model data) for the shown examples in Figure[T] In 
the figure, shadow areas represent the narrow emission lines, 
each double horizontal dashed lines represent the range of 
[/o — l,/o -I- 1], where /o = 0, 10, 20. ..80 representing the zero 
point for each spectrum are shown as solid horizontal line 
in the figure. The results shown in Figure [2] and in Figure [l] 
indicate that the elliptical disk model is better for double- 
peaked broad emission lines of 3C390.3, and further indicate 
that the probable existing hot spots and/or warped struc- 
tures have few effects on the measured flux density of broad 
Ha and few effects on our final results and conclusion, i.e., 
the hot spots and/or ward structures are weak for observed 
spectra of 3C390.3 around 1995, and elliptical accretion disk 
model is efficient and sufficient for 3C390.3 around 1995. 

Figure [3] shows the distributions of the seven disk pa- 
rameters, and the distributions of the parameter of log(x^) 
(the value of the summed squared residuals divided by the 
degree of freedom for the returned model parameter val- 
ues, which is one parameter as the residual to determine 
whether the theoretical model is preferred) for the 39 high 
quality spectra in upper three panels and for all the 66 
spectra in the other six panels. The final accepted disk pa- 
rameters are: inner radius ro = 216 ± 26-Rg, outer radius 
Ti = 1263 ± 70Rg, eccentricity e = 0.13 ± 0.04, inclination 
angle i — 29.57 ± 1.67°, orientation angle tjio = —23 ± 7°, 
emissivity power slope q = 1.65±0.26 (/(r) (x r^') and local 
broadening velocity a = 742±145km/s. Actually, we should 
note that some of the theoretical disk parameters, ro, ri, q, 
00, o", should depend on the strength of continuum emis- 
sion, in other words, the parameters are time dependent. 
However, there need to be long-term progressive changes in 
the ionizing fiux to change these variables. Moreover, from 
the results shown in Figure |3] we can find that there are 
tiny variations for these parameters, which indicates that 
the mean values for the theoretical disk parameters can be 
used to trace the actual physical disk parameters for 3C390.3 
around 1995. 

Before proceeding further, we compare our measured 
disk parameters with the previous results reported in the 
literature. In the sample of double-peaked emitters from ra- 
dio galaxies, Eracleous & Halpern (1994) gave the results 
about disk parameters for BLRs through circular accre- 
tion disk model (Chen et al. 1989) for double-peaked broad 
Ha observed in 1988 for 3C390.3 ( simple results can also 
be found in Figure 6 and in Section 4.5 in Sambruna et 
al. 2009), which are some different from our results deter- 
mined through elliptical accretion disk model (Eracleous et 
al. 1995). We know that in accretion disk model, full width 
at zero intensity (FWZI) of double-peaked broad emission 
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lines sensitively depends on inner radius of disk-like BLRs, 
smaller inner radius leads to broader FWZI. From the fitted 
results shown in Figure 4 in Eracluous & Halpern (f994), we 
can find that in order to better fit the broad wings of broad 
Ha, the inner radius should be some smaller than 380-Rg 
(the one listed in Eracluous & Halpern f994). Thus our inner 
radius ro = 2f6-RG should be reasonable. Peak separation 
sensitively depends on the outer radius, thus there are simi- 
lar results about outer radius in Eracluous & Halpern (f994) 
and in our paper. In order to better fit the cusp feature 
around blue peak, the elliptical disk model should be more 
efficient than the totally symmetric circular disk model, i.e., 
the parameters of eccentricity and orientation angle should 
be active. The separation between peak intensity and zero 
intensity of blue (or red) part of Ha should depends on the 
emissivity power slope. The larger distance between our in- 
ner radius and outer radius indicates that our emissivity 
slope should be smaller than 3 (the one listed in Eracluous 
& Halpern f994). Thus our emissivity slope q = f .65 should 
be reasonable. Line width of double-peaked broad line sensi- 
tively depends on inclination angle, thus it is clear there are 
similar value about inclination angle i ~ 30° in our paper 
and in Eracluous & Halpern (f994). Because the elliptical 
accretion disk model is more efficient for 3C390.3 than the 
circular accretion disk model, the effects from local broad- 
ening velocity should be lower. Thus our local broadening 
velocity a ~ 750km/s smaller than f 900km/s in Eracluous 
& Halpern (f994) should be reasonable. 

Besides the disk parameters, it is interesting to check 
the variations of line profiles of double-peaked broad Ha. As 
shown in Veilleux & Zheng (f99f ), there are significant vari- 
ations of relative flux ratio of the two peaks within the pe- 
riod from f974 to f988. Certainly, the variations can not be 
successfully explained by the elliptical accretion disk model, 
because the disk precession period of 3C390.3 based on the 
elliptical accretion disk model is about several hundreds of 
years as we should discuss in the following. However, during 
the period (about SOOdays) around 1995 included in AGN 
WATCH project, the variations of flux ratio of the two peaks 
are tiny as shown in Figure[4l In the figure, two kinds of fiux 
ratios are shown, the one with mean value of 1.01±0.02 (0.02 
is standard deviation for the mean value) is the ratio of the 
blue part to red part of broad Ha divided by 6564. 6lA (the 
theoretical center wavelength of Ha), the other one with 
mean value of 1.46 ± 0.06 (0.06 is standard deviation) is 
the intensity ratio of the blue peak to red peak. Similar re- 
sults for 3390.3 around 1995 can also be found in Dietrich 
et al. (1998). The more recent results about the flux ratio of 
the two peaks of 3C390.3 from 1995 to 2007 can be found in 
Shapovalova et al. (2010), there are tiny variations of the fiux 
ratio of the two peaks in the 5years period as shown in Figure 
13 in Shapovalova et al. (2010). Due to the tiny variations, 
we believe that there are tiny variations for disk parameters 
obtained from the observed spectra during 1995, and that is 
the reason why we select the spectra observed around 1995. 
In other words, the structures of BLRs of 3C390.3 are stable 
around 1995, which re-confirm that there are tiny variations 
of disk parameters shown in Figure [3] 

Last but not least, we should note the collected 66 
spectra in optical band for 3C390.3 are the observed spec- 
tra before the intercaliberation method (Van Groningen & 
Wanders 1992) applied to consider the effects from differ- 



ent observational instruments in different configurations as 
discussed in Dietrich et al. (1998). In this paper, we mainly 
focus on the disk parameters, rather than the flux densities 
of broad Ha which should be directly collected from AGN 
WATCH project with contamination being corrected. The 
main effect from intercaliberation method on the corrected 
line profile is the broadening velocity in the method. How- 
ever, the broadening velocity is only about tens of kilometers 
per second when intercaliberation method is applied, which 
is much smaller than line width of broad Ha. Thus, the ef- 
fects from intercaliberation method on the measured disk 
parameters through elliptical accretion disk model can be 
totally ignored. 



Before to finish the section, we consider the following 
question whether there are other values for disk parameters 
in the disk model applied to well fit observed double-peaked 
broad Ha. In other words, the question is whether the solu- 
tions to the disk parameters are unique. If the answer is yes, 
there should be one set of values for the disk parameters in 
the following mathematical procedure, and we will try to de- 
termine whether our mathematical procedure can be applied 
to determine the available parameter space. If the answer is 
no, the following mathematical procedure should be simple 
and succinct. To give one precise mathematical solution to 
the question is very difficult. We consider the question as 
follows. The observed double-peaked broad Ha is re-fitted. 
When procedure starts to fit line profile, one of the disk 
parameters is fixed to one value much different from the ac- 
cepted value above (half of the accepted value for the param- 
eter), and the other disk parameters with the same starting 
values in fitting procedure are free. Then the Levenberg- 
Marquardt least-squares minimization technique is applied 
to find the best fitted results for observed double-peaked 
broad Ha and find the best solutions for the disk parame- 
ters. Actually, the software package MPFIT in Markwardt 
IDL Library (http://www.physics.wisc.edu/~craigm/idl/ ) 
is used to perform the least-squares fitting and to find the 
best solution. In order to clearly compare the fitted results 
for different disk parameters, the value fluxfiti-fit2 is cal- 
culated, where fit2 means the best fitted results for accepted 
disk parameters shown in Figure [3] and fitl means the best 
fitted results for one fixed disk parameter with half of the 
accepted value for the parameter. Figure [5] shows the fit- 
ted results, and the corresponding values of fluxfiti-fit2- 
The apparent and large difference (residual to some extent) 
between fitl and fit2 indicates that solutions to disk pa- 
rameters based on accretion disk model are unique to some 
extent. Thus, in the following mathematical procedure to de- 
termine geometrical structure of emission line regions, there 
are fixed disk parameters. 



Based on the results above, through the disk param- 
eters, one kind disk-like geometrical structures of BLRs of 
3C390.3 can be well structured. Then it is very interest- 
ing to check whether the characters of observed light-curves 
can be reproduced through the theoretical disk-like BLRs of 
3C390.3, which is the main objective of the following section. 
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3 RESULTS ABOUT REVERBERATION 
MAPPING TECHNIQUE UNDER 
ELLIPTICAL ACCRETION DISK MODEL 

In the section above, one kind geometric structure of BLRs 
of 3C390.3 have been supposed by the theoretical elhptical 
accretion disk model (Eracleous et al. 1995). In this section, 
we will mainly check whether the supposed disk-like BLRs 
are valid to reproduce properties of observed continuum and 
broad line variability, through reverberation mapping tech- 
nique. 

In order to clearly build one physical geometrical struc- 
ture of disk-like BLRs, the length of one gravitational ra- 
dius {Rg) in physical unit should be firstly determined, 
i.e., it is necessary to determine the masses of central black 
hole of 3C390.3. Here the BH masses of 3C390.3, Mbh = 
4.68 X IO^Mq, is accepted. For 3C390.3, the two kinds of BH 
masses, both virial BH masses (Peterson et al. 2004) and 
BH masses from M-sigma relation (Tremaine et al. 2002, 
Ferrarese & Merritt 2001, Gebhardt et al., 2000, Giiltekin 
et al. 2009, Woo et al. 2010, Lewis & Eracleous 2006) can be 
found in literature, which are similar. Thus we do not worry 
about the accuracy of BH masses of 3C390.3, and accept 
Rg ~ 0.02631ight — days based on the BH masses. Through 
the disk parameters above, it is not difficult to build the 
physical geometry of the BLRs of 3C390.3 supposed by the- 
oretical model. Then it is interesting to check the results 
based on the reverberation mapping technique through the 
supposed disk-like BLRs as follows. 

First and Foremost, we accept the assumptions listed 
in Peterson (1993) for reverberation mapping technique, 1): 
the continuum emission is from one central source which 
is much smaller than BLRs. The assumption is efficient for 
3C390.3, although double-peaked broad line of 3C390.3 is 
assumed from accretion disk. The reported size of BLRs 
of 3C390.3, the time lag between observed optical balmer 
emission and observed UV X-ray emission, is about 20 ± 8 
light-days (Dietrich et al. 1998, Bentz et al. 2009), i.e., 
Rblrs ~ 400 — IOOORg, through cross correlation func- 
tion method. Based on the theoretical disk parameters for 
3C390.3, the flux weighted mean size of BLRs is about 
567Rg (simply calculated by (ro + ri)/(l. + q), q ^ 1.65 
is the line emissivity slope). The result indicates the UV 
emission region is much near to central black hole, which 
clearly can be treated as one point source. 2); both contin- 
uum emission and line emission are freely and isotopically 
propagating in the central volume, which can be confirmed 
by the small covering factor estimated by the appearance of 
no-absorption features around broad Ho. 3): the line emis- 
sions are in rapid response to ionizing continuum, which 
can be confirmed by the strong correlation between line lu- 
minosity and continuum luminosity found by Greene & Ho 
(2005) and by the much shorter recombination time (about 
one hundred seconds for standard BLRs) of balmer emission 
lines and much longer dynamic time (about several years for 
standard BLRs) than the light travel time across BLRs, as 
discussed in Peterson (1993). 

Besides, we structure the supposed theoretical struc- 
tures of BLRs of 3C390.3 as follows. More than 10^ tiny 
clouds (or so-called test particles, if we treat each tiny cloud 
with sufficient tiny size) are set in the N elliptical rings (the 
BLRs), each tiny cloud has its position (radius r and orien- 



tation angle (j)) and line intensity (/). The position in one 
elliptical ring for one tiny cloudy is created by. 



G [0, 2 X tt] 

_ r* X (1 + e) 
1 — e X cos((^) 



(1) 



where »"o ^ r* ^ ri and e are the pericenter distance and 
eccentricity of the disk-like BLRs of 3C390.3, r is the dis- 
tance from the tiny cloud to the central black hole. From 
inner boundary ro to outer boundary r\, the pericenter 
distance for the elliptical rings is evenly separated into N 
bins, log(ro) < log(r<.,i(* ~ l-.-N)) ^ log(ri) with step 
of (log(r-i) — log(r-o))/(A'^ — 1). Because most part of line 
emission is from inner area of BLRs, thus the parameter 
of logarithm of pericenter distance (log(r-)) is evenly sep- 
arated, rather than the parameter of pericenter distance 
(r), which provides more areas in inner part of BLRs. After 
that, the orientation angle is evenly separated into M bins, 
^ 0j(i = 1 . . . M) s; 2 X TT with step of 2 x 7r/(M - 1). 
Based on the r*,; and 4>ji the supposed theoretical disk-like 
BLRs can be separated into (A*' — 1) x [M — 1) tiny ar- 
eas {Ai.j{i — 1 . . . N,j — 1 . . . M)). As long as N is larger 
enough, the ionizing photo propagating time through each 
tiny area can be ignored, and each tiny cloud can be sim- 
ply treated as one point source. The line intensity from each 
area fij {i — 1 . . . N,j = 1 . . . A-I) can be directly calculated 
by elliptical accretion disk model, and the total line intensity 
F is the sum of ftj, 



fi,j ~ I I H{model)drd(j) oc P{ri,(f>j) 



(2) 



EE/- = ^ 



where H{model) is the integrand function used for accretion 
disk model as shown in Eracleous et al. (1995), g* is different 
from the emissivity slope in elliptical accretion disk model 
(however if the radius of each bin is uniformly created, then 
g* ^ 1.0 — g as shown in Figure O, P{ri,(j>j) is the func- 
tion depending on radius and orientation angle to calculate 
the line intensity of each tiny cloud in the ith elliptical ring 
which can also be found in Figure [SI Fi and F mean the line 
intensity from the ith elliptical ring and the total line inten- 
sity from total area. Here, we select A'^ = 60 and M = 400. 
A*' = 60 confirms that the light travel time from the i bin to 
i + 1 bin is much less than 1 day (especially for the bins in in- 
ner part of BLRs) , the standard date separation of our input 
light curve of continuum emission. Furthermore, we should 
note that in our procedure, the listed disk parameters above 
are taken as fixed values, based on the following consider- 
ations. On the one hand, there are tiny variations of line 
profiles of double-peaked Ha observed around 1995. Based 
on the observed Ha selected rom AGN WATCH, it is very 
difficult to find one reliable correlation between line profile 
variability and disk parameters (especially orientation angle 
0) through theoretical disk model. In other words, there is 
no information about the function of <?I)(t) where t is date and 
time. On the other hand, the fixed disk parameters through 
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tiny varied line profiles around 1995 lead to much simple 
mathematical procedure. 

Figure [6] shows the procedures above to build the struc- 
tures of BLRs of 3C390.3 based on the disk parameters, in- 
cluding the properties of equation (2). Once the tiny cloud 
meet the ionizing photo from central source, the line inten- 
sity of the tiny cloud is changed immediately as. 



conij(t — ij 



(3) 



where t is the date with uniform separation of 1 day as dis- 
cussed below, a is slope of the the corresponding correlation 
between luminosity of Ha and continuum luminosity, fij 
and coriij mean the line intensity from area A{i,j) and ar- 
riving continuum emission for area A{i,j). For AGN, a is 
the about 1 for QSOs in Greene & Ho (2005), and for low 
luminosity AGN in Zhang et al. (2008). However, a ~ 1 is 
not preferred for 3C390.3, Figure [7] shows the correlation 
between continuum luminosity at 5177A and luminosity of 
broad Ha for the spectra observed around 1995. The val- 
ues of luminosity of broad Ha and continuum luminosity 
at 5177A are collected from AGN WATCH, the contami- 
nation has been corrected. The spearman rank correlation 



coefficient is about 0.66 with Pnuii 



W 



Through the 



Levenberg-Marquardt least-squares minimization technique, 
the best weighted fitted result for the correlation can be cal- 
culated, 



8.92 X 10*^(- 



con(5177A) )0.26±0.03 ^ (^^ 

1044erg/s ' 6/ V ; 

Thus, we accept a ~ 0.26, not ~ 1. Actually, the different 
index of line-continuum luminosity correlation is due to the 
contamination of narrow lines around Ha. As shown in Di- 
etrich et al. (1998), the flux density of broad Ha includes 
contributions from narrow lines around Ha, and flux den- 
sities of narrow emission lines are constant for 3C390.3. If 
we define factor k as flux ratio of pure broad Ha to narrow 
lines, we should find that points with lower luminosity of 
Ha in Figure [7] have smaller values of k, points with larger 
luminosity of Ha in Figure [7| have larger values of k. If the 
contributions of narrow emission lines are corrected, the in- 
dex for the correlation shown in Figure [7| should be steeper 
than 0.26 and near to 1. 

Figure[8]shows one simple output line curve of broad Ha 
based on one test input light curve of continuum emission 
described by one delta function through the disk-like BLRs 
determined elliptical accretion disk model for 3C390.3. Be- 
cause of the extended BLRs, the delta function is delayed 
and expanded. Light traveling from the central black hole 
(we assume that the ionizing photos are from central point) 
to the innerest boundary of the BLRs is about 5days, from 
the inner boundary (~ 2167?g) to the point with longest 
distance (~ 1640-Rg, larger than the pericenter distance 
1263-Rg for the out boundary, because the BLRs are ellipti- 
cal disk-like) from central black hole is about 39 days. Thus 
in the figure, it is clear that 5 days after the burst of contin- 
uum, the line intensity starts to vary, and then 39 days after 
the burst of continuum emission, the variation of line inten- 
sity ends. Because the inner part of BLRs emits most of the 
line intensity of broad Ha, thus the line intensity of total 
Ha is changed from strong to weak with passage of time as 
shown in the figure. Furthermore, it should be interesting 



to check the effects of different geometrical structure on the 
response output broad Ha, although we have accepted there 
are fixed and stable disk parameters in section above. We 
accepted the disk parameters for 3C390.3 under the stochas- 
tically perturbed accretion disk model in Flohic & Eracleous 
(2008, one circular disk part with bright spots), ro = 450-Rg, 
ri = UOORg, a = 1300km/s, i = 27° and <? = 3, in spite of 
the not well fitted results for line profiles with these param- 
eters (not good fitted results for broad wings). Under the 
model, the output Ha is also shown in Figure [S] It is clear 
that the response to input continuum described delta func- 
tion is much different for elliptical disk structure and for cir- 
cular disk structure. The different response for different geo- 
metrical structures indicates the different structures of BLRs 
of AGN should be discriminated through the light curves of 
continuum emission and broad emission lines emission, espe- 
cially when there are HOMOGENEOUS and COMPLETE 
information for the light curves. Actually, based on the ob- 
served light curves of 3C390.3 selected from AGN WATCH 
project, there is no way to discriminate different geometrical 
structures of BLRs as what we have shown in the following 
results. The results shown in the figure indicate that our 
procedures above DO work, and effects of extended disk-like 
BLRs on light curve of emission line can be clearly described 
by our procedures. 

Last but not least, the observed light curve of contin- 
uum emission at ~ 5177A collected from AGN WATCH is 
used as the input light curve of continuum emission. The de- 
tailed description about the light curve can be found in Di- 
etrich et al. (1998). There are 70 data points included in the 
light curve from MJD2449517 to MJD2450068. We first cre- 
ate one new light curve of continuum emission with date sep- 
aration of Iday by linear interpolation method applied to ob- 
served light curve of continuum. Then we check the expected 
output light curve of broad Ha under the supposed elliptical 
disk-like BLRs. The results are shown in bottom-left panel 
of Figure [S] It is clear that the expected output light curve 
of broad Ha is consistent with the observed light curve of 
broad Ha, the linear correlation coefficient for the correla- 
tion between output light curve of Ha and observed light 
curve of Ha is about 0.85 with Pnuti ~ 0. Then through the 
input light curve of continuum and the output corresponding 
light curve of broad Ha, it is interesting to estimate the size 
of BLRs by so-called cross correlation function (CCF). Here, 
the common interpolated cross-correlation function (ICCF) 
(GaskeU & Sparke 1986, GaskeU & Peterson 1987, Peterson 
1993) is applied to quantify time lag between continuum 
emission and broad lines emission. We do not consider the 
z-transfer discrete correlation function (ZDCF, Alexander 
1997, Edelson & Krohk 1988, White & Peterson 1994) any 
more, because results from ZDCF are excellent agreement 
with the results from ICCF (Peterson et al. 1991, 1992, 2004, 
Kaspi et al. 2000, Bentz et al. 2010). Furthermore, the cor- 
responding results based on the disk parameters shown in 
Flohic & Eracleous (2008) are also calculated, and shown in 
Figure [9l It is clear that although the example shown in Fig- 
ure[S]indicates the information from light curves can be used 
to confirm which structure should be more preferred, there 
is no way to find enough evidence to determine which struc- 
ture (elliptical disk-like BLRs or circular disk-like BLRs) is 
preferred, due to the incomplete and in- homogeneous infor- 
mation provided from the observed light curves. Thus, we 
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mainly compare the calculated time lags between continuum 
emission and broad line emission, under the two structures. 

Figure[9]shows the CCF results from input light curve of 
continuum and the output light curve of broad Ha, the dis- 
tribution of the size of BLRs through bootstrap method (a 
common Monte-Carlo method to estimate the uncertainty of 
parameter. Press et al. 1992, Peterson et al. 1998). Detailed 
descriptions about CCF method to estimate size of BLRs 
(time lag between line emission and continuum emission) 
and bootstrap method to estimate uncertainties of measured 
time lag can be found in Peterson (1993). From the results 
shown in the figure, we can find that the estimated size of 
BLRs through the input light curve of continuum emission 
and output light curve of broad Ha under elliptical accre- 
tion disk model is also well consistent with the one esti- 
mated by the observed light curves of continuum emission 
and broad line emission, Rblrs ~ 16.8 ± 2.11ight — days. 
We should note that the size of BLRs is the distance from 
broad Ha emission line region to optical continuum emission 
region, thus the value is some smaller than ~ 201ight — days 
estimated by the CCF function from light curve of broad 
Ha and light curve of x-ray/UV emission shown in Diet- 
rich et al. (1998), Bentz et al. (2009) etc., because there 
are probable 5days lag between optical continuum emis- 
sion and x-ray emission (Dietrich et al. 1998). Certainly, 
we also calculate the results based on the circular disk 
model shown in Flohic & Eracleous (2008), which should 
be some different from the results for elliptical accretion 
disk Rblrs ~ 18.4 ± 2.11ight — days. However, there is no 
enough evidence to against the circular disk-like BLRs for 
3C390.3, besides the best fitted results for observed line pro- 
files. If there should homogeneous and complete observed 
light curves, the clear final decision should be given. 

In order to further confirm that the disk-like BLRs de- 
termined by elliptical accretion disk model is reliable for 
double-peaked emitter 3C390.3, we also check the CCF re- 
sults for disk-like BLRs with different disk parameters, in 
spite of the best fitted results for observed double-peaked 
broad balmer line. Figure [10] shows the effects of disk pa- 
rameters on the measured size of BLRs through input con- 
tinuum emission and output Ha emission based on the el- 
liptical accretion disk model as what we have done above. 
It is clear that inner radius and emissivity power are the 
two main parameters which have apparent effects on the 
theoretically measured size of BLRs of 3C390.3. In the fig- 
ure, we only show the effects from inner radius, emissivity 
power, eccentricity and inclination angle of disk-like BLRs, 
because there are few effects from the other disk param- 
eters. Certainly, more reliable evidence to determine and 
confirm model parameters should depend on future more 
detailed and homogeneous observed light curves of contin- 
uum emission and broad line emission. The results indicate 
that the measured disk parameters for BLRs of 3C390.3 are 
reliable to some extent, and different geometrical structures 
of BLRs based on different disk parameters should lead to 
different theoretically measured size of BLRs from size of 
BLRs through observational results. 

The results above indicate that the elliptical accretion 
disk model can be applied to best fit the observed spectra, 
i.e., the model can provide fine velocity structures for BLRs 
of 3C390.3. Furthermore, the elliptical accretion disk model 
can provide fine geometrical structures which can be applied 



to reproduce characters of observed light curve of broad Ha. 
Thus elliptical accretion disk model which can be applied to 
fit characters of line profiles and variations of broad double- 
peaked Ha is appropriate to double-peaked emitter 3C390.3. 



4 DISCUSSION AND CONCLUSION 

As one well-known double-peaked emitter, based on the disk 
parameters measured above, to check the expected line pro- 
file of double-peaked broad balmer emission line of 3C390.3 
in different years/periods is very interesting. Although the 
more recent observed line profiles of double-peaked broad 
Ha of 3C390.3 in Jan. 2007 can be found in Shapovalova et 
al. (2010), the observed spectrum in Sep. 2000 by HST STIS 
is collected, because the reduced spectrum in 2000 is public 
and conveniently collected from the website of Multimission 
Archive at STScI (MAST) (http://archive.stsci.edu/) (HST 
Proposal 8700, PI: Prof. Andre Martel in Space Telescope 
Science Institute). The detailed description about observa- 
tional technique for the spectrum around 2000 by HST can 
be found in Popovic (2003). Here, we mainly check whether 
the observed line profile of broad balmer emission lines in 
2000 can be reproduced through theoretical elliptical ac- 
cretion disk model based on the disk parameters measured 
from spectra observed around 1995 for double-peaked emit- 
ter 3C390.3. 

Before proceeding further, the precession period should 
be determined. Through the obtained disk parameters, the 
relativistic precession period of the elliptical disk-like BLRs 
around central black hole of 3C390.3 can be simply calcu- 
lated (Weinberg 1972), 



2 X TT 



2 X TT 



Gtt X G X Mb h 



''Ail- 



') 



70years (inner radius) 
370years (outer radius) 



(5) 



where 'A' is the semi-major axis length and e is the ec- 
centricity. Although, the precession period for outer ellipti- 
cal rings is very long, the short precession period for inner 
elliptical rings should lead to some apparent variations of 
double-peaked broad emission lines from 1995 to 2000. 

We simply assume that orientation angle is the unique 
parameter to vary from 1995 to 2000 in the elliptical ac- 
cretion disk model. As what we have done in section 3, 
the elliptical disk-like BLRs are separated into 60 rings. It 
is clear that orientation angle of each ring in 2000 should 
be some different due to different semi-major axis length 



for each ring, 



± 2 X TT X 



where ± 



means the rotating direction of disk-like BLRs is clock- 
wise or anticlockwise for observer. After the orientation an- 
gles of the 60 rings around 2000 are determined, line pro- 
file expected in 2000 can be re-constructed by the sum 
of 60 expected double-peaked components from 60 rings, 
as what we have done in Section 3. Here the theoretical 
line profile in 2000 is mainly based on the mean spec- 
trum around 1995 in AGN WATCH project. The mean ob- 
served spectrum shown as thin solid line in Figure [TT] is 
created by PCA (Principal Components Analysis, so-called 
Karhunen-Loeve Transform) method applied for all the 66 
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observed spectra of broad Ha around 1995. The conve- 
nient and public IDL PCA program 'pca_solve.pro' writ- 
ten by D. Schlegel in Princeton University is used, which 
is included in the SDSS software package of IDLSPEC2D 



| http://spectro.princeton.edu/| . Commonly, the first prin- 
cipal component represents the mean spectrum. From the 
mean spectrum, it is clear that there are similar flux ratio 
(~ 1.46) of blue peak to red peak of broad double-peaked 
Ha and similar flux ratio (~ 1) of blue part to red part with 
those mean values shown in Figure U 



Figure [TT] shows the probable observed line profile 
of broad double-peaked Ha of 3C390.3 around 2000, and 
shows the comparison between the mean observed spectrum 
around 1995 and observed spectrum around 2000. Based 
on the theoretical accretion disk model, through the results 
shown in top-left panel in Figure [TTl we can find that due to 
precession period, after about 5 years, the flux ratio of red 
peak to blue peak should be changed. In the figure, fiux den- 
sity at blue peaks of all the line profiles have been normalized 
to 1. If rotating direction of disk- like BLRs is clockwise to 
observer, the expected fiux ratio (1.36) of red peak to blue 
peak of spectrum around 2000 should be smaller than the 
one (1.46) of observed spectrum around 1995. If rotating di- 
rection of disk-like BLRs is anticlockwise to observer, the 
expected fiux ratio (1.68) in 2000 should be larger than the 
one (1.46) in 1995. Through the observed spectrum in 2000, 
the flux ratio of blue peak to red peak is about 1.79 ± 0.18. 
Thus, if anticlockwise rotating disk-like BLRs are applied for 
3C390.3, the variations in observed spectra in different years 
can be naturally explained, although the applied model is 
over-simplified. The some larger value 1.79 than theoreti- 
cal expected 1.68 are perhaps due to the not so accurate 
precession period, the probable existence of varying of hot 
spots and/or warped structures, different local broadening 
velocity (which should lead the part around red peak to be 
some flat), etc.. The not so consistent red peak positions 
for spectrum around 2000 and mean spectrum around 1995 
is probable due to the variation of emissivity power slope 
during the 5 years from 1995 to 2000. 



Finally, a simple summary is listed as follows. Based on 
the observed spectra around 1995 for 3C390.3, the elliptical 
accretion disk model is applied to fit the observed double- 
peaked broad Ha collected from AGN WATCH project. 
Then based on the disk parameters, the formed geomet- 
rical structures of BLRs of 3C390.3 are applied to check 
the reverberation mapping results for 3C390.3. The disk- 
like BLRs supposed by theoretical elliptical accretion disk 
model can well re-produce the observational results about 
reverberation mapping technique. Thus the disk-like BLRs 
are preferred for 3C390.3. Furthermore, we check the effects 
of precession of elliptical disk-like BLRs on the observed line 
profile. After about 5 years, the line profile of broad double- 
peaked Ha observed in 2000 by HST can be well expected by 
the theoretical elliptical accretion disk model, based on the 
disk parameters measured through spectra observed around 
1995, which confirms the elliptical accretion disk model is 
appropriate to double-peaked emitter 3C390.3 . 
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Figure 1. Best fitted results for double-peaked broad Ha in rest wavelength by elliptical aceretion disk model. Thin solid line represents 
the observed spectrum, thick solid line represents the best fitted results for double-peaked broad Ho. The vertical dashed line shows the 
position for unexpected absorption features around 6500A for two examples marked with '49860ce' and '49984ce'. The observed MJD 
date of each spectrum is shown in the right side of the figure. Shadow areas represent the ranges for narrow emission lines. 
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Figure 2. The residuals (yobs ~ Vmodeh observed data minus expected model data) for the shown examples in Figure [T] Shadow 
areas represent the ranges for narrow emission lines, each double horizontal dashed lines represent the range of [/o — l,/o + 1]) where 
/o = 0, 10, 20... 80 representing the zero point for each spectrum are shown as solid horizontal lines in the figure. 
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Figure 3. Distributions of the disk parameters of disk-like BLRs for 3C390.3. The top three subfigures show the distributions of 
eccentricity, inclination angle and the parameter of log(x^) for the 39 spectra with high quality and without unexpected absorption 
features around 6500A. Then the other six subfigures show the distributions of the other parameters for all the 66 spectra, inner radius 
ro, outer radius ri, emissivity slope q, orientation angle </>o, local broadening velocity a and the parameter of log(x^)- The mean value 
of each parameter we accepted is shown in each subfigure. 
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Figure 4. The variations of the flux ratio of blue peak to red peak (shown as open circles), and the flux ratio of blue part to red part 
of broad Ha (shown as solid circles). The dotted lines represents the mean values of the two kinds of flux ratios. 



3C390.3: disk-like BLRs near central BH 13 




6700 5800 



Figure 5. Best fitted results for double-peaked broad Ha with one fixed parameter with half of accepted value for the disk parameter. 
In the top panel, thin solid line in histogram mode represents the observed line profile marked with '49870dr' in AGN WATCH project, 
the thick solid line shows the best fitted results with disk parameters having values shown in Figure [S] thin solid line represents the best 
fitted results with disk parameter ro = 127.1-Rg, thin dotted line represents the best fitted results with disk parameter r\ = 670.8-Rg, 
thick dotted line represents the best fitted results with disk parameter sin(i) = 0.25, thin dashed line represents the best fitted results 
with disk parameter q = 0.90, thick dashed line represents the best fitted results with disk parameter a = 364fcm/s, thin dot-dashed 
line represents the best fitted results with disk parameter e = 0.052, thick dot-dashed line represents the best fitted results with disk 
parameter <f>0 = —13 in degree. The other seven panels show the corresponding values of fluxfiti_fit2) where fit2 represents the best 
fitted results with disk parameters shown in Figure |3] 
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Figure 6. Top left panel shows the toy structures of the elliptical disk-like BLRs of 3C390.3. The solid line represents the elliptical ring 
with pericenter distance rg ^ r^ ^ ri, the innerest one with pericenter distance of tq, the outerest one with ri. The central black hole is 
located at the origin of the coordinates (one focus point of the elliptical ring). Dotted line represents the sphere surface of ionizing photos 
at one time. Because the BLRs are elliptical disk-like, thus at one time, the ionizing photos arriving in BLRs DO NOT affect all the tiny 
clouds located in one elliptical ring. Top right panel show the toy showing of the line intensity of broad Ha from different bins of radius. 
Solid line represents the observed Ha, dotted line represents the line intensity from ~ 216Rg (the inner boundary) to ~ 480i?G, dashed 
line represents the line intensity from ^ 480-Rg' to ~ 740Rq, dot-dashed line represents the line intensity from ~ TAORq to ^ lOOO-Rc 
double-dot-dashed line represents the line intensity from ~ IOOORq to 1263Rg (the outer boundary). The vertical dotted line represents 
the center wavelength of Ha, 6564. 6lA. Bottom left panel shows the correlation between line intensities from different bins {Fi{r)) and 
the corresponding pericenter distances of the bins. Here the bins are created uniformly, thus F{r) oc r'^'^"'" oc r^-0-&5) ^ Bottom right 
panel shows the properties of line intensities (fij) of the 23541 data points in sixty elliptical rings with pericenter distances from ro to 
ri. The bottom two panels show the properties of Equation (2). 
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Figure 7. The correlation between luminosity of broad Ho and continuum luminosity at 5177A for 66 observed spectra around 1995. 
Solid line represents the best fitted result, Lh^ gc (L "j.)"'^®^. The dotted line is the one found by Greene & Ho(2005). 
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Figure 8. The response output broad Ha based on the input continuum emission described by one delta function under the elliptical 
disk-like BLRs and circular disk-like BLRs for 3C390.3. The top line shows the continuum emission. The bottom solid line shows the 
corresponding response output broad Ha under elliptical disk-like BLRs, and dotted line shows the results under circular disk-like BLRs. 
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Figure 9. Top left panel shows the light curve of continuum, open circles are the observed data points selected from AGN WATCH, 
solid line represents the input light curve of continuum with separation of Iday. Bottom left panel shows the light curve of broad Ha, 
open circles are observed values, solid line represents the output light curve under the elliptical disk-like BLRs, dotted line represents the 
output light curve under the circular disk model shown in Flohic & Eracleous (2008). Top right panel shows the CCF function (maximum 
coefficient about 0.99) for observed light-curve of continuum emission and output light curve of broad Ha. Solid line represents the result 
for elliptical accretion disk model (peak value around 16 days), dotted line shows the result for circular disk model in Flohic & Eracleous 
(2008) (peak value around 18 days). Bottom right shows the distribution of time lag between observed continuum emission and simulated 
broad Ha emission by bootstrap method. Solid line is for elliptical disk model, dotted line is for circular disk model in Flohic & Eracleous 
(2008). 
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Figure 10. The effects of disk parameters on the theoretically measured size of BLRs. Top-left panel shows the effects of inner radius, 
top-right panel shows the effects of emissivity power, bottom-left panel shows the effects of eccentricity, bottom-right panel shows the 
effects of inclination angle. 
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Figure 11. To compare the observed line profile by HST in 2000 and mean observed line profile around 1995. Thin dotted line represents 
the mean spectrum around 1995 and thick solid line represents the observed line profile by HST in 2000. Two vertical dot-dashed 
lines mark the positions of red peak and blue peak of double-peaked broad Ho. Top left panel shows the theoretical results. Solid line 
represents the mean observed double-peaked line profile around 1995, thick dashed line represents the expected line profile around 2000 
with clockwise rotating disk-like BLRs, thick dotted line represents the expected line profile around 2000 with anticlockwise rotating 
disk-like BLRs. 



